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AERODYNAMIC FORCES ON A VIBRATING 
UNSTAGGERED CASCADE* 

By H. Sohngen 


The unsteady aerodynamic forces, [based on two -dimensional incom- 
pressible flow considerations] , are determined for an unstaggered cascade, 
the blades of which Eire vibrating in phase in an approach flow parallel 
to the blades. 


INTRODUCTION 


In the theory of axial turbomachines, the aerodynamic forces acting 
on a vibrating cascade are of interest for investigations of vibrations. 
These forces were determined in reference 1 for the • case where the 
spacing of the cascade Is small. Information regarding an arbitrary 
staggering does not yet exist; we shall investigate here the simplest 
case of this kind, namely an unstaggered cascade with straight profiles, 
the blades of which vibrate in phase. We do not presuppose to limit 
our consideration to pure bending or torsional vibrations; rather, arbi- 
trary periodic deformations of the blade are admitted and their lift 
distribution is determined. This includes also the case where the flow 
approaching the blades is unsteady. We do assume, however, that the 
blades are not under static load, that the fluid is incompressible and 
frictionless, and that the amplitudes of vibration are small. 


LIFT DISTRIBUTION AND LIFT 


Assume the approach flow toward the cascade to be parallel to the 
blades with the velocity V. Let the velocity with which the blade 

point x' moves in the y' direction be given by g(x') x e iv ^. This 
velocity is the same for all blades . We visualize the blades of the 

cascade as covered by vortices of the density y(x' ) . It must be 
noted, in addition, that free vortices separate from the trailing edges 

of the individual blades in proportion as the total circulation about 

— — — — — ■ — — — 
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the blade in question varies with the time, although with the inverse 
sign. This is stipulated by the theorem of the conservation of circu- 
lation in the total space. In a linear theory, for the amplitudes of 
vibration to which we here limit ourselves, we may assume that these 
free vortices move away with the basic flow. If we denote accordingly 
their density by 


V 



(x* > 2 / 2 ) 


there applies at the trailing edge (x 1 = 2 / 2 ) the relationship 
d P 1 / 2 ( ivt a « „ iv(t-2/2V) 

~ / 7 (x’)e dx* + V 7 0 e = 0 (l) 

dtJ_ z / 2 0 

by which the density of the free vortices is connected with the varia- 
tion of the circulation about the blade. 

The vortices situated on the blades 7 (x’) and the free vortices 
7 q have to be determined - with consideration of the relationship (l) - 

in such a manner that the velocity component in the y* direction, 
induced by all the vortices together at the blade point x', is equal 
to that with which the blade point moves in this direction, that is, 
eq ual to g(x ' ) . A row of vortices, which lie on the blades at the 
point x' =0 and each of which has the circulation 7 , induce a field 
whose y' component on the blades at the point x' is equal to 


OO ^ f 

ct, a 

e - e 

This expression is easily obtained by superimposing the fields of the 
individual vortices and by taking into consideration since this is an 
axial-flow problem, that there can be no velocity induced far upstream of 
the cascade, at development. This makes it necessary to superimpose 
another transverse velocity on the summed row of vortices. Using the 
expression given above, we then obtain as a boundary condition the inte- 
gral equation 



NACA TH 1412 


3 




x'-£’) 


& 


J(x’-r) -§( 


x’-f) 


Zfi r e 

a -'l/2 


Vt i 






£(x'-t') -J(x’-I’) 

e a - e a 


<34’ = g(x') 


for which we may also write 


n l/2 

f 7(ft ') 

a J-l / 2 


^Lx 1 

- OCr 


fr x ' ^ 

e - e 


d& 1 - 



a 





2jt x t 2itt t 

a "aT 5 

e - e 


d£ 1 = g(x’) 


( 2 ) 


From this equation and equation (l), we have to determine the vortex 
density 7(x’)- 


Rowe ver , it is not so much the vortex density which is physically 
interesting but the lift distribution ^£p(x'), that is, the pressure 
difference between positive pressure side and suction side. We conclude 
from the unsteady-pressure equation that between lift distribution and 
vortex density there exists the relationship 

X * 

P(x’) = r(x’) + i £ f 7(l 1 )d| 1 (3) 

V J -1/2 

where 

Z®(x’) = pVP(x’ 

Thus, the lift distribution may be calculated as soon as the vortex 
distribution is known. 


In determining the lift distribution for the single blade, one 
usually starts from the integral equation which is satisfied by the 
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lift distribution. Therefore, we shall indicate the latter for the 
present case. It reads 


o 1/2 

4 P(|')K(x' - |')cL6' - g(x\) (|x'|<l/ 2 ) (4) 

v -Z/2 


■ Z/2 


with the kernel 


K(x') = -± 


£x’ 


a 2.x * 
a 


- e 


-£x’ 

ct 


+ i 


If 

VaJo 




g(x'-u’) 


g(x'-u') - g(x'-u') 

e a - e a 


du' 


(5) 


If we note that, because of (5), 


r(x') = p(x') - i % [ P(S') 

-1/2 


-i^(x'-S') 


Aft' 


we can show that the integral equation (4) is identical with (2) if 7 q 
satisfies equation (l) . 

However, since the kernel of the integral equation (4.) can he less 
easily seen through than that of (2), it appears more suitable to deter- 
mine the lift distribution not directly by inversion of equation (4) 
but first to solve equation (2) with respect to 7 (x'), with an as yet 
undetermined constant 7 q. Thus we obtain 7 (x') as a function of g(x’) 

and 7 q. If we introduce it into equation (l), a linear equation for 7^ 

results from which 7^ may be calculated as a function of g(x' ) . 

Thereby 7 (x') also is determined, and the lift distribution can be 
calculated. 

For solving the integral equation (2), one's first impulse is to 
try to reduce it back to the integral equation 


*/.: 


f(i) 


d| = G(x) 


( 6 ) 


known from lifting-surface theory. For this purpose, we set 


■ 2x' 

1 


X 


(7) 
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and 


*1 = 



TX 


•] 


(8) 


with 


t = Til /a, s = Sinh r c = Cosh r 

Furthermore, we stipulate that the primed coordinates shadl always be 
dimensional quantities; those without primes, in contrast, dimensionless 
as they result by means of equation (7) - Coordinates with the sub- 
script 1 shall always result from the coordinates without primes by 
means of the substitution (8) . In a manner which cannot be misinter- 
preted, we write 


g(x') = g(x) = gCx-jJ 

and correspondingly also, other functions . If f(x) is a given func- 
tion, we always understand by f j_(x^) the expression 


*l(xi) = 


f(*L) 


SXq + C 

Thus, the integral equation (2) then is transformed into 


n r 1 (S-,) 

St J x~ - ~ ii d| i = ~ 7 o w i (x i ) “ g i (x i } 


where w^(x]_) denotes the function described by 


( 9 ) 


Wl(xi) - as/j. 


.-toi d| 

*1 - 61 


(10) 


and 


cd = vlf 2V 

is the so-called reduced frequency. Furthermore, the relationship (l) 
is equivalent to 


„ . ffis 1m 

7q = -i _ e 


J 7i(x 1 )dx 1 


T 


(H) 
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Prom the integral equation (9), there follows, (reference 2) 


?1 ( X1 ) - — s /lZjS Z g- ifc) * . g ^ 6j ) Ml (12) 

jr ^1 - x x 2 J -1 y 1 - X X 2 x l " &i 


where the constant C at first may still he chosen arbitrarily. The 
condition of a smooth outflow at the trailing edge requires that 


0 ■ 'fl yrr§ [ 2r ° Wi N + 

Hence 7^(xj_) then takes the form 

’•<■>> - 1 pl/I 0 ; 

If we finally take into consideration that 

C = T y-^x^dx-L 
w -1 

there results for the constant y^ the value 

To gl<Xl) \/rrt 


?0 - 


B(oi, 


where 


b( “’ t) - T? - 2 


(13) 


(14) 


(15) 


(16) 


(17) 


is a function which depends only on the reduced frequency and the geo- 
metrical dimensions of the cascade. 

Thus the problem is solved in principle. The function g(x') is 
given; with it, the constant 7 q may be calcinated from equation (l6), 

equation ( 14) then yields the vortex distribution 7-j_( x l) > and the lift 
distribution results from 
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P(x) = r 1 (x 1 )(sx 1 + c) + iui | J 7i(&3)fl4i ( l8 ) 

It is now only a matter of writing this expression in a form more suit- 
able for numerical calculation. It will be seen, in particular, that, 
whereas the influence of the function - that is the influence 

of the free vortices . - enters into the vortex distribution in a rela- 
tively complicated manner, its effect on the lift distribution is of a 
considerably simpler structure . 


THE FUNCTION T(o, t) 


Before deriving a simpler representation of the lift distribution, 
we shall give a series expansion for the function B(o},t). From (10) 
and ( 17 ) , it follows that 

, \ e~^ 1 r 1 /I + x, r m e-* 1 * 

B( “’ T) " d| 3X1 

If we interchange the integration sequence and take into account that 




1 I 1 + x i 


1 - x L x x - I-l 


= -it 


*1 +1 
4l - l 


- 1 


ih > i ) 


there results 


-itu r * 0 
B(co,t) = — — + / 

ito J 1 




i 


Si + 1 


il - 1 


- li 


d| 


( 19 ) 


Furthermore it is convenient to split the function B(ca, t) into two 
parts . We let 


and 



( 20 ) 


Bi(a>, t) 



( 21 ) 
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Then 


B(cd,t) =||B 0 (oi,T) + B 1 (cD,T)j’ (22) 

Particularly simple representations exist for these two functions if 
a Jl = °o, that is, r = 0, so that we have the case of the single blade. 
Since, = it we see immediately that 

B 0 (co,0) = -IEq^Co)) 


where B^ 2 ) is the Hankel function. Furthermore, we can easily show 
that 


B^O) * -H^^Cco) 


We shall now derive representations for these two functions which 
are convenient especially for moderate values of a/l. We have 


B 0 (oi, t) 



If we understand by C the path of integration represented in figure 2, % 

the integral taken over it vanishes. Since the expression under the 
square root sign takes equal values on the two parallel paths, there 
follows, with 


ca' 



(l - e 27KD ')Bo(av 


■) = 


i 2 ae -to-r 

3T 


fo 


2a / 1 




\j (e. lT7 ^ - l)(e iTTl - e“ 2T ) 


dT! 


= i 


2 s e “isi-T r 2 * 
Jt T J 0 


m'-a 


- l)t“ - e- 2T ) 


dd 




n=0 
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However (ref. 5, p. 138 ) 


L 


( 05 ' -in) 






i-a 


and therefore. 


ri 2lt 

Jo 


-ii + ' -in-i^} d 

L ______ dd 


^ 2 sin d/2 

rfn + ^ + iio'^r^-n - 5ai') 


Bo(ai,T) = 


a»‘ 


JgL - 


-2t) -im 
e /e 


t( 1 + 21to f )r(-n - ici' )r^n + ^ + 1m ’j Sinh iOi) 


n=0 


c^e 


-2nr 


where 




E ( v ~ * k 0 ') 

V=1 v + ^ + im 1 j 


a Q = 1 


Because of 




< 1, this expansion converges so well for the spacing 
factors customary in turbomachines <x/l < 1, that is, t 2; it, that we 
generally can get by with only the two first terms. The function B]_(<d,t) 
may be calculated in exactly the same manner, and we obtain 


B^m, t) = - 


*?!*.■** 


iT*(l - ln)')r^i+ w\sinhioD 




-2nr 


Therefore, for the function which, later on, will be the only one 
of interest 


I(b,t) = E j--' ^ 

B 1 + B o 


( 23 ) 


which, for t = 0, is transformed into the function known from the theory 
of the single wing (ref. ^) , namely, 
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T(*o) = 


H^ 2 ^(o)) - 1Hq( 2 )(*o) 
H-^^oo) + IBq^ 2 ^ (to) 


there follows the representation 


T(< 


\ . 23*0' /, -2t\ 

.00 , t) =1 — 1 1 - e ) 

1 + 23*0' 


oo 

? 


a n e 


-2nT 


• rrfi? £(■ • »)-•“' 


(24) 


Prom this we see that for spacing factors a/l < 1 with an error of 
not quite 1 percent we may set 


1 - T(co, t) = 


2i*u ' 


1 •+ 23*o ' 


(24') 


LIFT DISTRIBUTION AND LIFT 


If we write the vortex distribution in the form 


7i(*i) 


2 

rt 



7q v i(xi) + g 1 (x 1 ) 
' X l 2 


cLx-l + 


2 

it 




V^l) + »l(*l) a , 
( X 1 - h )^ 1 - t * 


there follows (ref. 5, Al) 




(25) 


(26) 
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with 


where 


0 <; arc cos x., <; it 


w l(jj»l)^l = J' ®l(^l)^l (27) 

and C has the significance given above; namely 

-ico 

< 28 ) 

In order to calculate the function WgCx ^) , we start from the 
function 


w i u i> o) - sbl 


.-®6 dt 


X 1 - *1 


x^ < 1 


For this function, w i( x i> ia>) = w i( x i) - If we first assume that Rcr > 0, 
there applies 

/ e '° S to (Sl - ^l) 4 * 

- *1 
- -Us^ Ml - * x ) 4- 4 r e -° 8 — 3- — — ail 

2as a x - x-j_ d| 

7, „-cr T P“ -erg sfl-, - x n ) + SX-i + c 

= — ln(l - x-,) + — / e 5 -All U i d|> 

2 as a ea |i - x^ 


l e “ ff 


Hence it follows that 


r> 1 

/ wp(u3_, cr) du^ = ’ 7 ~~ e ^(1 - x i) - (sx-,_ + c) ^ w-lCx-^c) 


If we let cr -» ira, there follows 


W 2^ x l) = 2 ^ ^ 111(1 ‘ Xl) ‘ ^ (SX1 + C) V I (xi) 
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If we further consider the integral relationship 


/; T^\-P , -.x - Al - i arc cos J 

J - 1 ( x i - Si) vA - il 5 L * y 


with 


0 < arc cos it 

we obtain from (l8) for the lift distribution the representation 


* I 


g 2 (vi)dvi 


(xi - V X )^1 - V;L 2 


r op ax i + c) \/i-j^ - 8 i/ rT ^ 5 lf z 1 
L v 1 j J -1 


dvi 


in which only the integral over the function w-j_ remains to be calcu- 
lated. For it there results 

2 r 1 ^i(^i)a ... = _L_r 1 __i r“ e-^ 

11 - 1 x/ 1 - 1 51(16 v/ -1 '\jl - V-Zj 1 V x - dS dVl 


=-— r 
1 


-ln)| 


- 7 o <% 

v /i 7 7ri 


I n 

2as 


®o^ T ) 


and we obtain, 'because of (l6) 
- 1 w 1 (v 1 ) 


‘ 7 ° dn = Kl ' T(1 “' T)] gl(Xl) \/r I S 


dxi 
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where T is the expression defined "by (25) . There follows finally the 
ultimate representation for the lift distribution 


p(x) - 1 (sxj. * JsM d 6l ♦ 

^ 2 \ji - x, 2 f at, + 

T ‘ v 1 (-1 - H) yrn? 1 

| E 1 - T(cd,t)] |( 8X1 + c) - 8 ^ L ' x l 2 | 


(29) 


which for t -» 0 is transformed into the representation known for the 
single wing (refs. 4, 5)- For the lift we obtain from ( 29 ) 


= r 2 4, ax' = 1 pvi a r 1 - p(x) . ax x x e ivt 
J-7/2 2 T J_l8X 1 +C 


■ PVl f -! 




1 g C3 (x 1 )(l - C - S Xl ) 


(c + SX' 


i ) V rr 


e 


ivt 


( 30 ) 


BENDING VIBRATIONS 


As an example of application, we shall consider the case where the 
blades of the ring perform bending vibrations. If A is the amplitude 

of the deflection, 1 oj 2V ^ e^ v "^ is the velocity with which each blade 
point moves in the direction of the y* axis. Therefore, 

g(x’) = io}2V y = 6 = const 

L 


and hence there follows from the definition given above : 


g 1 (x 1 ) = 


and g 2 (x^) = — ln(s Xl + c) 
s 


SX-]_ + c 



NACA OM l4l2 


l4 


If we enter these into (30) , we obtain for the lift the expression 


A = «pV 2 A 



T) 


— 1 — + as 2 — In Cosh 1 1 (31) 

1 + e T J 1 + e T ts 2 


which for t -» 0 is transformed into 


*PV 2 A ’joj 2 - ioj(l + T(cjo) Jj” 


That is the value known for the single wing (ref. 4) . If'wq assume, 
on the other hand that r is large, and' limit ourselves to a linear 
theory in a/l, we obtain, in agreement with reference 1 


A = pV 2 A | - 2 ioJ 


(32) 


If we use for spacings a,/l <L 1 the approximation (24’) for 1 - T, 
the lift may be written in the form 

A = pV 2 A y jto 2 fi(u>',T) - 2 ia>f 2 U',T)j* (33) 

with 


f 1 (^p 1 , t) 

f 2 (a>’,T) 


= — In Cosh ^ 
r 2 


= e T - e~ T 
e T + 1 


e T - e" T 


r(l + W 2 )(l + e T ) 2 

W 2 1 

(l + W 2 )(l + e T )J 


where these two functions assume the value 1 for a linear theory in a/l. 
Hence we can see that for spacings a/l < 0.5; that is, for t values 
^ 2it, the linear theory yields the left term which is proportional to 
the velocity; thus the damping, with an error of not quite 1 percent. 

In contrast, for the term which is proportional to the deflection, that 
is, behaves like a spring force, conditions are considerably more unfa- 
vorable. This term is very inadequately included by the linear theory. 
We must note, however, that in vibration calculations this term is 
opposed by a mechanical spring force which is generally very large com- 
pared to the corresponding term of an aerodynamic force so that the 
deficiency has only a slight effect. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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Figure 1. - Cascade parameters. 




Figure 2. - Path of integration C. 
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